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We use the quasiclassical theory of superconductivity to calculate the electronic contribution to 
the thermal conductivity. The theory is formulated for low temperatures when heat transport is 
limited by electron scattering from random defects and for superconductors with nodes in the order 
parameter. We show that certain eigenvalues of the thermal conductivity tensor are universal at low 
temperature, A;_bT ^ 7, where 7 is the bandwidth of impurity bound states in the superconducting 
phase. The components of the electrical and thermal conductivity also obey a Wiedemann-Franz 
law with the Lorenz ratio, L{T) = n/aT, given by the Sommerfeld value of Ls = {■k^ I ?>){k b I eY 
for A;_bT ^ 7. For intermediate temperatures the Lorenz ratio deviates significantly from Ls, and 
is strongly dependent on the scattering cross section, and qualitatively different for resonant vs. 
nonresonant scattering. We include comparisons with other theoretical calculations and the thermal 
conductivity data for the high Tc cuprate and heavy fermion superconductors. 
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I. INTRODUCTION 

In a normal metal at sufficiently low temperatures 
the electrical and thermal conductivities are determined 
by the scattering of electrons by lattice defects. The 
electrical conductivity approaches a constant, while the 
heat conductivity, k{T) ~ T, is related to the electri- 
cal conductivity by Sommerfeld's result for the Lorenz 
ratio, k/((tT) Ls = {ir^ /3){kB/ef . In fact the 
Wiedemann-Franz (WF) law is frequently used to esti- 
mate the phonon contribution to k by subtracting off the 
expected electronic contribution, aLsT, from the mea- 
sured heat conductivity.^ 

Superconductivity has dramatic effects on the electri- 
cal and thermal conductivities in conventional (s-wave) 
superconductors.^^® In particular, the WF law is violated 
by the formation of a coherent ground state and the open- 
ing of a gap in the excitation spectrum everywhere on the 
Fermi surface. In this article we investigate the behav- 
ior of the heat current for superconductors with an order 
parameter of reduced symmetry for which there are gap- 
less excitations even for the pure superconductor. Such 
superconducting states have been argued to exist both in 
the cuprates^ and the heavy fermion systems. In par- 
ticular, it is widely believed that the order parameter for 
the heavy fermion superconductor UPts vanishes on a 
line in the basal plane, pf;, = [e.g., see Ref. 11]. Sim- 
ilarly, one of the leading candidates for the cuprates is 
the Big or dj;2_y2 state, a singlet order parameter with 
lines of zeros at the Fermi surface positions, pf^ = ipfy 
[for a review see Ref. 9]. 

For a clean superconductor with an order parameter 
that vanishes along a line on the Fermi surface the den- 



sity of states is linear in the excitation energy, A(e) ~ 
Nf e/Ao for e < Aq."'^^'"'^^ However, this spectrum is al- 
tered by a random distribution of impurities. ^'^'^^ A new 
energy scale, 7, develops below which the density of 
states is approximately constant and non-zero at zero en- 
ergy. The energy scale 7 is interpreted as the bandwidth 
of quasiparticle states bound to impurities. These 
impurity bound states develop below the superconduct- 
ing transition and are coherent superpositions of particle 
and hole excitations. Such states are formed by the con- 
structive interference of particle- and hole-like excitations 
that undergo Andreev scattering from the variations of 
the order parameter that occur as a result of unconven- 
tional pairing and potential scattering by the impurity. 
For an order parameter with a line of nodes, the band- 
width, 7, and the density of Andreev bound states at zero 
energy, A^(0), are finite for any finite concentration of im- 
purities, riifnp 7^ 0.^*^ The energy scale, 7, and the density 
of states, A^(0), depend on both the impurity concentra- 
tion, riimp, and the scattering phase shift, 60- Thus, 7 
provides a crossover energy scale, below which the trans- 
port properties of an unconventional superconductor are 
dominated by the Andreev bound states. For excitation 
energies above 7 the transport properties are determined 
primarily by the scattering of continuum excitations. 

The electrical conductivity for a superconductor with 
an order parameter that vanishes along a line of nodes 
was shown theoretically by Lee^"'^ to have a universal 
limiting value, (Tq = lira^^o a(u),T = 0) ~ e^NfV^r^, 
where ta — 7i/(7rAo) is a universal transport time that 
is independent of either the concentration or the scat- 
tering phase shift. This result was obtained for a two- 
dimensional (2D) dj;2_y2 order parameter of the form 
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A(pf) = Ao cos(2i^) and an isotropic 2D Fermi surface, 
p-f = (cos 1^, sin 1^). Calculations show that the univer- 
sal value of the conductivity is due to the phase space for 
"optical" transitions within the band of impurity-induced 
Andreev bound states. ^^'^^ The phase space requirements 
for a universal limit are determined by the variation of 
the order parameter near the nodes; A(i9) ~ — i9node 
for a line node in 3D or a point node in 2D, and A(i9) ~ 
— i^node)^ for a point node in 3D. The crossover to the 
universality limit occurs for ksT ^ 7 and kuj ^ 7. Thus, 
the universal limit is most easily realized in the strong 
scattering (unitarity) limit. Since there is considerable 
evidence that some of the heavy fermion superconduc- 
tors have an order parameter with a line of zeroes, the 
experimental confirmation of the universal result for the 
conductivity would provide an important test of the ar- 
gument (based on the Kondo lattice model) that impu- 
rity scattering in heavy fermion metals is in the unitar- 
ity limit. Unitarity scattering by impurities has also 
been invoked in the high Tg superconductors in order to 
reconcile the often observed dependence of the pen- 
etration depth at low temperatures with a rf-wave order 
parameter. 

Because electromagnetic fields penetrate only a dis- 
tance of the order of the London length into the super- 
conductor, surface effects can complicate the determi- 
nation of the bulk conductivity. On the other hand, 
the heat current is unscreened and provides more di- 
rect access to the bulk excitation spectrum. In this 
article, we investigate the low temperature behavior of 
the thermal conductivity tensor for unconventional su- 
perconductors with line and point nodes in the order 
parameter. One of the issues we address is whether or 
not the universal behavior of the electrical conductiv- 
ity extends to the electronic heat conductivity at low 
temperature. We show that the components of the elec- 
tronic thermal conductivity tensor k corresponding to 
quasiparticles in the vicinity of the line nodes are de- 
termined by the same scattering rate as the electrical 
conductivity and are universal in the limit T ^ 0. Fur- 
thermore, the WF law is obeyed for the ratio of the uni- 
versal electrical and thermal conductivities in the limit 
ksT 7. However, a significant temperature depen- 
dence of the Lorenz ratio occurs over the temperature 
range, T < Tg, even with purely elastic scattering. The 
universal values for both the electrical and thermal con- 
ductivity result from the cancelation between two factors: 
(i) the density of Andreev bound states, which is pro- 
portional to 7, and (ii) the reduction of phase space for 
scattering of gapless excitations, which is proportional to 
7~^, leading to an estimate for the thermal conductivity, 
K ~ Nf{j/Ao)klTvf{n/j) ~ NfvfklT{n/ Ao), which 
is independent of the defect density or scattering phase 
shift (again for n~^p > ^q). Perhaps the most surpris- 
ing result is that the ratio of the universal values for the 
thermal and electrical conductivity gives the Sommer- 
feld value for the Lorenz ratio, k/uT 2^ Ls = ^(ks /e)^ ■ 
Thus, the differences in the coherence factors that deter- 
mine the conductivity tensors, k and a , do not affect the 
Lorenz ratio L(T) = n/aT for ksT ^ 7 and kuj ^ 7. 

For temperatures above the crossover energy, ksT ^ 7, 



the Lorenz ratio, L(T), deviates significantly from the 
Sommerfeld value. Furthermore, we find that the tem- 
perature dependence of the Lorenz ratio is very sensitive 
to the scattering phase shift. For nearly resonant scatter- 
ing [60 ~ ''■/2] L(T) is larger than Ls , except in a narrow 
region near Tg associated with a coherence peak in the 
electrical conductivity for very clean superconductors. In 
the opposite limit of weak scattering the Lorenz ratio is 
less than Ls except for the (exponentially small) region 
ksT ^ 7. Thus, measurements of L(T) might be useful 
in distinguishing weak and strong impurity scattering. 

The rest of this paper is organized as follows. In Sec. 
II we derive an equation for the thermal conductivity of 
a Fermi liquid; the formulation includes unconventional 
pairing and the effects of scattering by a random dis- 
tribution of defects. In Sec. Ill we evaluate the ther- 
mal conductivity tensor in the limit T ^ for several 
models of the order parameter: (i) the even-parity Bi^ 
state in tetragonal symmetry, i.e., the dj;2_y2 state with 
A(pf) = Ao(p^^ ~ P]yy^ ^^'^ hexagonal systems, (ii) 
the even-parity Ei^ ("hybrid-I")^^ state with A{pf ) = 
2Aopj^{pjx + ipfy), (iii) the odd-parity Ai^ ("polar") 
state with A{pj) = Aozp/^, and (iv) the odd-parity 

E2„ ("hybrid-II") state A{pi) = ^Aozpf,{pf^+ipfyf. 
The dj;2_y2 state has been discussed extensively as a 
model for the high Tg cuprates, and the latter three or- 
der parameters have been discussed as models for the 
low-temperature superconducting phase of UPts. For 
the odd-parity states (cases iii and iv) z specifies the 
quantization direction for the spins; the pairs are in the 
triplet spin state, [ti + it); relative to the z direction. 
The formulation and many of the results are applica- 
ble to more general forms of anisotropic pairing with ze- 
roes in the order parameter. A few other models are 
also discussed. In Sec. IV we discuss the leading or- 
der finite temperature corrections to the thermal and 
electrical conductivities. For an order parameter with 
a line of nodes, or a point node in which the gap opens 
quadratically, the thermal conductivity tensor has com- 
ponents that are universal in the limit T ^ 0, and of 
the form k = LsaT{l + 0[T'^/j'^]), exhibiting both the 
WF law for T ^ and the deviations that develop for 
ksT ~ 7 ^ Aq. In Sec. V we present numerical re- 
sults for the thermal conductivity and the Lorenz ratio 
over the full temperature range below Tg, and compare in 
Sec. VI the results with low-temperature measurements 
of the thermal conductivity for several cuprate and heavy 
fermion superconductors. 



II. QUASICLASSICAL TRANSPORT 
COEFFICIENTS 

The microscopic theory of superconductivity was de- 
veloped by Bardeen, Cooper & Schrieffer^*^ at about the 
same time that Landau published the microscopic basis of 
his transport theory of normal Fermi liquids. These two 
theories were combined into what was called by Larkin 
& Ovchinnikov^^ the quasiclassical theory of supercon- 
ducting Fermi liquids. The quasiclassical theory of su- 
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perconductivity is a complete theory of the low-energy 
properties of fermions in the superconducting state. It 
was developed by Eilenberger,^^ Larkin & Ovchinnikov,^^ 
and Eliashberg,^'* and it covers essentially all thermody- 
namic, electrodynamic, transport and collective proper- 
ties of superconductors. In this section we start from 
the basic equations of the quasiclassical theory of super- 
conductivity, and derive the electrical conductivity and 
the electronic contribution to the heat conductivity of 
anisotropic superconductors with unconventional pairing 
and scattering by random defects. 

The central physical objects of the quasiclassical the- 
ory are the quasiclassical propagators, which obey quasi- 
classical transport-like equations. We give a brief inter- 
pretation of their physical meaning, and establish our 
notation. We use Keldysh's formulation of nonequi- 
librium Green's function theory,^® and introduce three 
types of propagators: advanced (A), retarded (R), and 
Keldysh (K). These three quasiclassical propagators are 
4 X 4-matrices whose components describe the quantum- 
mechanical internal degrees of freedom of electrons and 
holes: the spin and particle-hole degrees of freedom. Par- 
ticle and hole excitations are incoherent in the normal 
state, whereas the superconducting state is characterized 
by quantum coherence between particles and holes, which 
is the origin of persistent currents and other nonclassi- 
cal superconducting effects. The quasiclassical propaga- 
tors describe the quantum statistical state of the internal 
degrees of freedom. Nonvanishing off-diagonal elements 
in the particle-hole index indicate superconductivity. A 
standard notation for the matrix structure of the propa- 
gators (and the self-energies) is^® 

x_( a^'+g''-^ (f" +P'-^)i'ry\ 
^ ~ \ ■ ffX ^ fX A X -X ^ 

with X E {i?, j4. A'}. The 16 matrix elements of are 
written in terms of four spin scalars {g-^ , g^ , , / ) 

and four spin vectors {g-^ , g^ , , ) . All matrix 
elements are functions of the Fermi momentum pj , the 
position R, the excitation energy e, measured from the 
chemical potential, and the time t. The diagonal spin 
scalars g^ , g^ contain the spectral and statistical in- 
formation for spin-independent quantities like the heat 
current density, while the spin vectors, g^ , carry the in- 
formation on the spin magnetization, spin currents, etc. 
The off-diagonal terms and characterize the su- 
perconducting state; a finite value of indicates singlet 
pairing, a nonvanishing / implies triplet pairing. 

The redundant information provided in the definition 
of the quasiclassical propagators in Eq. (1) can be elimi- 
nated with the very general symmetries^® 

g^{p};e) = g^{ps-e)\ g^{pj-e) = g_^{-ps--e), (2) 

a'^iP! ; e) = a^^ip! ; 0* > a^^iP! ; = a^i-P! ; -0 > (3) 
/^(p/;0 = /^(p/;0% /^(P/;0 = /^(-p/;-0, (4) 
tiP! ; = P'iP! ; 0* , P'iP! ; = -f^i-P! ; -0 • (5) 

The electrical current density is obtained from the 
scalar part of the Keldysh propagator, the Fermi velocity. 



vj , and the density of states per spin, Nj , at the Fermi 
level, 

ie(i?,t) = 2#; j dpj j ^,[evj{pj)]g^{pj,R;e,t), (6) 

where J dpj{. . .) stands for a normalized integral over the 
Fermi surface. Similarly, the heat current density has the 
form, 

UR,i) = 2Nj j dps j ^. [evsmg^'ipj^ii-e^i). (7) 

For weak disturbances from equilibrium the current re- 
sponse is linear in the applied field. In this paper we are 
interested in the low-frequency, dissipative part of the 
electrical current response, defined by the conductivity 
tensor, 

8]e=^ -E^, (8) 

where (T = lin\^^o Re (t(w, T),^^ and the thermal con- 
ductivity tensor, defined by the linear response to a small 
temperature gradient 

8l = -K.VT. (9) 

In the rest of this section we develop the linear response 
equations for *k and (T from the quasiclassical theory. The 
analysis and notation closely follows that for the current 
response to an EM field given in Refs. 38 and 22. The 
advanced, retarded, and Keldysh propagators are calcu- 
lated from quasiclassical transport equations, 

[eh - <Te.t - a^'^ , g""'^] „ + ivj ■ Vg^'^ = , (10) 

and 

(era - agxt - o-^) ° 9^ - 9^ ° {eh - &ext - o"^) 

-a^ o'g"^ +'g^oa^ +ivf -Vg^ = (11) 

where all propagators and self-energies depend on the 
Fermi momentum pj , the position R, the excitation en- 
ergy e, and the time t. We use the units k = ks = I, 
unless explicitly stated. The o-product stands for the 
following operation in the energy-time variables 

(ao'h^ (pf,R;e,t) = 

ei(9?9?-9?9,^)a(^^ ^ R. i)6(^^ ^ R. ^ (12) 

and the commutator [a, &]o is defined by aoh — hoa. The 
transport equations are supplemented by the normaliza- 
tion conditions, ^^'^^ 

g^'"^ o g^'^ = -TT^ l, (13) 
g^ o gK + gK o g^ = . (14) 

The quasiclassical transport equations (lO)-(ll) together 
with the normalization conditions (13)-(14) and the equa- 
tions specifying the self-energies, (T^, are the fundamen- 
tal equations of the Fermi liquid theory of superconduc- 
tivity. They are the generalization of the Boltzmann- 
Landau transport equation to the superconducting state. 
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The transformations and approximations used to derive 
transport equations are based on a systematic expansion 
to leading order in the small parameters of Fermi liq- 
uid theory, e.g., ksTc/Ef, Tilou/Ej, Ti/tEj, etc. The 
accuracy and predictive power of Fermi liquid theory is 
intimately connected with the smallness of these param- 
eters. 

The quasiclassical self-energy terms, (5"^, in the trans- 
port equations describe interactions between quasipar- 
ticles with phonons, with impurities, and quasiparti- 
cles with each other. We consider the low temperature 
transport properties of superconductors with unconven- 
tional pairing under conditions where inelastic scattering 
by phonons and quasiparticles is negligible compared to 
scattering from random defects. 

The quasiclassical self-energies depend on interaction 
vertices which are phenomenological parameters of the 
Fermi liquid theory of superconductivity. We consider 
the weak-coupling limit^® in which electronic pairing in- 
teractions are described by the vertices Vi - , and Vi ~ , , 

■' VsVj VsVj ' 

for the spin-singlet and spin-triplet interactions, respec- 
tively. The mean-field pairing self-energies are given by 

A«.^(p;,i?;t)= ldp/V^^^yf^(p/,R;e,t), (15) 



de 



A^'^(p;,i?;t)= jdp/Vi^pyf^ip/,R;c,t), (16) 

(17) 



= = . 



The effects of a random distribution of impurities are 
described in Fermi liquid theory by an electron-impurity 
vertex, u(pf,p^'), and the impurity concentration riimp-^^ 
The impurity self-energy is proportional to the single im- 
purity t matrix. 



t^(pf,Pf,R;e,t), (18) 



where t-^ are obtained by solving the integral equations. 



P'^(p; , p/, R; e, t) = u{pf , p/) l+Nf J dp/' u{pf , p/') 

g'^'^ip/', R; e, t) o i^'^{p/',p/, R; e, t) , (19) 

i''(pf,p/,R;e,t) = Nf j dp/' {""(pf , p/' , R; e,t) 

og'^ip/', R- e, t) o i^{p/',p/, R; e, t) . (20) 

The coupling of quasiparticles (with charge e) to an 
electric field is given by 



O-ext = -- Vf ■ Afa . 

c 



(21) 



terms of the deviations from local equilibrium specified by 
a thermal distribution function with a local temperature. 



$o(i?) = [l-2/(e;T(i?))] = tanh 



2T{R) 



(22) 



The local equilibrium Keldysh propagator and self- 
energy are determined by the retarded and advanced 
functions and the thermal distribution function. 



O-f = (T^ O $0 - $0 O &^ 



(23) 
(24) 



with retarded and advanced propagators that are given 
by the solutions of 



-.R,A -.R,A 



-.R,A 

9o 



0. 



-ttM 



(25) 
(26) 



and the self-consistency equations (15)-(20). Note, that 
the o-product reduces to matrix multiplication for the 
local equilibrium functions. The self-energy includes the 
mean-field order parameter, A(pf,R), and the impu- 
rity self-energy, a^^(pf , R;e), which has both diagonal 

(S^p) and off-diagonal (A^^) components in particle- 
hole space. 

In this paper we consider only superconducting states 
which are "unitary", i.e., the equilibrium mean-field or- 
der parameter satisfies 



A(pf,Rf 



|A(P/)|M 



(27) 



where |Ap stands for either the spin scalar product A A 
or the spin vector product A • A. The unitary condi- 
tion restricts us to even-parity, spin-singlet pairing or to 
odd-parity, spin-triplet states without spontaneous spin 
polarization. The odd-parity states considered in this 
paper are unitary states that do not break time-reversal 
symmetry in the spin degrees of freedom. However, time- 
reversal symmetry may still be broken by the orbital mo- 
tion of the Cooper pairs, which is the case for the Ei^ and 
E2ti ground states that we consider. 

The local equilibrium solutions to Eqs. (25)-(26) for 
unitary states in unconventional superconductors [with 
fdpf A(pf) = 0] are^o 



On (pf , R;e) = —TT — , , 

° ' ' 7|A|2 - (g«,^)2 

e^'^ipj , i?; e) = e - 1 Tr [htf^^^ips , R; e)] 



(28) 
(29) 



where A{q, to) is the vector potential describing the trans- 
verse field E = —A. In order to calculate the conductiv- 

c 

ity we must solve the transport equations to linear order 
in the perturbing field. We first linearize the transport 
equations in the perturbations of g^ and (T^ from their 
equilibrium values. In the case of the heat transport the 
perturbation is the temperature gradient VT. Thus, the 
linearization of the transport equations is carried out in 



These equilibrium functions are inputs to the linearized 
quasiclassical transport equations. The quasiclassical 
transport equations and normalization conditions are 
solved to linear order for the deviation of the propaga- 
tors from their local equilibrium values, Sg-^ (pf , R; e,t) = 

g'^ {pf, R', c,t) — g^{pf, R', e). The technical steps used to 
decouple the retarded, advanced and Keldysh functions, 
and for inverting the linearized transport equations are 
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outlined in the Appendix, and the solution for 6g is 
given in Eq. (A20). In the following we use the general 
solution for 6g^ and the self-consistency equations for 
the impurity self-energy and order parameter to obtain 
formulas for the electrical and thermal conductivities for 
a superconductor with an unconventional order parame- 
ter. 



III. ELECTRICAL AND THERMAL 
CONDUCTIVITIES 



Here we consider a superconductor with anisotropic 
singlet pairing or unitary triplet pairing, and discuss the 
electrical and thermal conductivities in the long wave- 
length limit q ^ and at T ^ 0. For simplicity we 
assume isotropic impurity scattering. In this case the 
first order corrections to the current response functions 
from the impurity self-energy, 6aimp, and the order pa- 
rameter, 6A, vanish for all listed pairing states, except 
for the polar state (ii) with current fiow along the c- 
axis}^'^^'^^ Self-energy corrections corresponding to the 
excitation of collective modes of the order parameter, 6 A, 
also vanish in the limit q ^ [cf. Refs. 43 and 44]. The 
self-energy corrections are the 'vertex corrections' in the 
language of the Green's function formulation of the Kubo 
response function.'*® If we can neglect vertex corrections, 
we obtain expressions for the electrical and thermal con- 
ductivities that depend only on the equilibrium propa- 
gators and self-energies, and the external perturbations. 
For spin singlet states the thermal conductivity obtained 
from Eqs. (7) and (A20) becomes 



Kij(T) 



47rT2 



de 



sech 



V2TJ 



7r2C|(p/;e)2-KL>«(e)2 

-Li{P};e)f^{pj;e) + ^^]. (30) 

In case of spin triplet pairing the off-diagonal spin scalar 
Green's functions in (30) have to be replaced by f^f^ 

f^-f^- Note, that only the anomalous part of the propa- 
gator (see Appendix) contributes to the thermal conduc- 
tivity. The retarded and advanced parts drop out after 
taking the trace and applying the normalization condi- 
tion, i.e., Tr (/^'^ V (/^'^ = 0. Physically, this means that 
the deviation of the quasiparticle distribution function 
due to a thermal gradient contributes to the heat cur- 
rent, whereas changes in the quasiparticle and Cooper 
pair spectrum do not. Equation (30), combined with the 
equilibrium propagators, impurity self-energy and order 
parameter, is the basic result for the electronic contri- 
bution to the thermal conductivity tensor. Note, that 
we have used the short-hand notation: C;|?^'^'"(f>f ; e), and 
D^'^^'^^e), for the functions in Eq. (30) at w = 0. It can 
be shown that Eq. (30) for Kij reduces to the same ex- 
pression for the thermal conductivity as reported previ- 
ously by Schmitt-Rink it et al.,^'* Hirschfeld it et al.,'*®''*^ 
and by Fledderjohann & Hirschfeld;'*^ except that these 
authors appear to have dropped the term from the 
impurity self-energy, which, however, vanishes in both 
Born and unitarity limits. 

Similarly, the electrical conductivity obtained from 
Eqs. (6) and (A20) for a spin singlet state is given by 



Re aij(ui,T) 



e^Nf 



2t:uj 



de 



tanh 



(-) 

\2TJ 



tanh ( — - I 
\2TJ 



Re- 



Cl{P}\e,uj) 



TT^Cl(pf;e,Loy + D<'_(e,Loy 



[g^ipf ■ e.)g^(pf ; e+) + f^(pf ; e.)f^(pf ; e+) + tt^) 



(31) 



where £± = e ± w/2. For triplet pairing simply re- 
place /^(e_)/o(e+) /o(f-) • /o(f+)- This result was 
obtained earlier for magnetic scattering in conventional 
superconductors,'*® and for electron-phonon and impurity 
scattering in strong coupling superconductors,^® and for 
the in-plane conductivity of layered superconductors. 
The formula for the conductivity reduces to the well- 
known result of Mattis and Bardeen for dirty, s-wave 
superconductors, 2 and to the result derived in Ref. 49 for 
the frequency and temperature dependence of the con- 
ductivity of a d-wave superconductor with lines of nodes 
in the order parameter. 

In deriving Eqs. (30) and (31) we have made use 
of several relations which are consequences of the gen- 



eral symmetries of the propagator (and self-energy) in 
Eqs. (2)-(5), the parity of the order parameter, and the 
specific symmetry of the equilibrium Green's function 
g^ = —(/if, which follows directly from Eq. (28). The ba- 
sic functions defining the self-energy and response func- 
tions obey the following symmetries'®'^® 

?^(e) = ?«(e)* , ?^(e) = -g«(-e), (32) 

C^ipj;c,i,) = Cfipj;c,i,r , (33) 

Clipf;e,cu) = Clipf;e,-cur , (34) 

L>^(e,w) = L>f(e,w)* , (35) 

£)«(e,w) = -£)«(£, -w)* . (36) 
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A. Wiedemann- Franz Law for T ^ 



In the limit T ^ and w ^ the occupation factors [tanh((e+^)/2T) — tanh((e— ^)/2T)] and sech^(e/2T) confine 
the e-integrals in Eqs. (30)-(31) to a small e region (of order T or uj). Assuming that there exists an energy scale, 
e* ^ T, on which the propagators and self-energies vary, we can set e = in the slowly varying parts of the integrands 
and obtain 



Nf . f 2 ( ^ \ 

Re (Tj', (w ^ 0, T ^ 0) = lim — — ^ sinh — — / rfe sech i — i 

U)^0 /TTW 2i J 



f ■ / 2 / 

— smh / rlf sech I 



\2TJ 



CliPf) 



7r2C|(p/)2 + L>« (0): 



[9^{Pf)9^{Pf)+k^{Pf)f^{Pf) + ^')], (37) 



and 



where the energies and frequencies are fixed to zero in the arguments of the propagators and self-energies. Using 
the symmetry relations (2)-(5) and (32)-(36) and eliminating the advanced and anomalous functions in Eqs. (37) and 
(38), we find 

Re^,,(c.^O,T^O) =^ Jdesech'[^) J dpf [vf.Vf,]^^^^ , (39) 



'^.•(^-0) = A/ c^ee2sech^(^) / dp^v^.v,:] ^f^y (40) 



We 



sd the normalization 



iditi' 



/f(P/)/o^(P/) 



9o(pfy 



-TT , in addition to the symmetry re- 



lations, in order to put the momentum integrals in Eqs. 
(39) and (40) in identical form. It is useful to write our 
final results in terms of an average Fermi velocity and a 
tensor, j , which incorporates all of the coherence effects 
of superconductivity at T ^ into an effective transport 
scattering time. The energy integrals are standard, so the 
conductivities for a system with D dimensions reduce to 



Re aij (w ^ 0, T ^ 0) = e^-Nf v] nj 



D 
3 D 



(41) 
(42) 



where vj = J dpj\vj{jtj)[^ , and the effective transport 
time is defined by the tensor 



D 



b/.(P/)^/,j(P/)]?''(0) 



Ho)2 



(43) 



For an isotropic normal metal one has ?^(0) = i/2TQ, 
where tq is the quasiparticle lifetime due to impurity 
scattering in the normal state. The transport lifetime 
in the normal state reduces to tq for isotropic impurity 



scattering, i.e., Tij = To6ij . Note, that Eq. (43) is ap- 
plicable to the normal state because the key assumption 
in deriving Eqs. (39)-(40) was that T is small compared 
with e* , where e* is the energy on which the propagators 
and self-energies vary. Thus for the normal state e* ~ £"/ , 
while for the superconducting state e* ~ 7, where 7 is the 
impurity bandwidth. In some respect the impurity band 
states form a new low temperature metallic state deep in 
the superconducting phase. This analogy is strengthened 
later in the paper when we calculate the temperature cor- 
rections to the transport coefficients using a Sommerfeld 
expansion. However, the "metallic" band of impurity 
states has other properties that differ significantly from 
those of conventional metals. The special features of the 
impurity induced metallic band refiect the reduced di- 
mensionality for the phase space of scattering and the 
energy dependence of the particle-hole coherence factors, 
which define the impurity induced band. These two fea- 
tures lead to (i) universality for T ^ transport coef- 
ficients for excitation gaps with line nodes or quadratic 
point nodes, and (ii) to the temperature dependence of 
the Lorenz ratio for elastic scattering at < T < T^. 

However for T ^ 0, we note that one obtains from Eqs. 
(41) and (42) the Wiedemann-Franz law with Sommer- 
feld's result for the Lorenz number of an unconventional 
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superconductor at very low temperatures, Kij = Ls <Tij T , 
with 




We emphasize that Eqs. (41), (42) hold for gapless su- 
perconductors in which the leading contribution to the 
transport current is that from quasiparticle excitations 
with energies T < j Tg. For superconductors with a 
gap the number of quasiparticle excitations at low tem- 
perature is activated, cx exp (— Aq/T), and the transport 
coefficients in the limit T Tg are not described by Eqs. 
(41), (42), and furthermore do not obey the Wiedemann- 
Franz law. 



isotropic impurities also obey the WF law. Although 
scattering by impurities is expected to be the dominant 
scattering mechanism at low temperatures, it is worth 
noting that these universal results also hold for electron- 
electron scattering or electron-phonon scattering, as long 
as vertex corrections can be neglected. 



IV. LOW TEMPERATURE LIMITS FOR 
SEVERAL UNCONVENTIONAL PAIRING 
STATES 

A. Zero temperature limit 



B. Universal Limits for d-wave pairing 

Lee has shown that for T ^ the electrical conduc- 
tivity of a d-wave superconductor is universal, ^"'^ 



(T||(w ^ 0,T^ 0) ~ e^NfV^TA, 
with the universal transport scattering time 



1 



TrAn 



(45) 



(46) 



which is independent of the scattering rate of quasipar- 
ticle excitations. 

Our quasiclassical calculation reproduces Lee's result 
for isotropic systems, pf = pf (cos ip, sim^), and the stan- 
dard model for rf-wave pairing, i.e., A(pf) = Aq cos(2i^). 
We obtain the universal result in the limit, |?^(0)| ^ Aq, 
where Aq is the maximum gap. The transport lifetime 
tensor from Eq. (43) for the standard rf-wave model re- 
duces to Tii = r 6ii , 



dip cos^(ip)j^ 



27r 



[A(^)2 + ^2]- 
where j is the width of the impurity band, 
(7(A2+7^)-i/2) 



7 



cot2,5o + (7(A2 -K72)-i/2y 



(47) 



(48) 



and < ... >= /o^" If (•••)) and r„ = nimp/(irNf). For 
a given impurity concentration this bandwidth is largest 
in the limit of unitarity scattering, e.g., for dj;2_y2 pair- 
ing, 7 ~ ^7rAort,/2, while in the Born limit, 7 ~ 
4Ao exp(— 7rAo/2rt,). In either case we have 7 ^ Aq, 
otherwise pair-breaking by impurity scattering effectively 
suppresses the superconducting transition. For 7 ^ Aq, 
Eq. (47) reduces to r ~ l/(7rAo), and we obtain Lee's 
universal result in Eq. (45). 

Since, according to Eq. (43), the same transport life- 
time determines the electronic contribution to the ther- 
mal conductivity, it too becomes universal in the limit 
T ^ 0, i.e., independent of the impurity scattering rates. 
The universal, low-temperature limits for the electri- 
cal and thermal conductivities of a standard rf-wave su- 
perconductor with A cx cos(2i^) and scattering due to 



We evaluate the eigenvalues of the k's in the T ^ 
limit for the specific pairing states listed in Table I for 
uniaxial superconductors. For heat fiow in the basal 
plane the two eigenvalues of Kij are identical, and we 
drop the subscripts for these directions and write K|| = 
Kaa = Kill,, and Kj_ = K^c for the c-axis transport coeffi- 
cient. 

First, consider the dj;2_y2 pairing state. For a cylin- 
drical Fermi surface and the heat flow in the a&-plane, 
the angular average in Eq. (47) for the transport time 
reduces to i • 4 • f , |^ 7^ (A^ 72)-3/2. There are 
four line nodes and the integral is reduced to one quad- 
rant containing one node. For 7 ^ Aq, the integral is 
dominated by the region near the node, thus, we can ap- 
proximate A(i^) ~ iJ,Aoip' , where ip' = 7r/4 — ip and the 
parameter /j, measures the slope of the gap at the node, 
i.e., /J, = (l/Ao)|(iA(i^)/(ii^|(^^^^^. This slope parameter 
also determines the low-energy density of states [cf. Ref. 
50]. The result for the low-temperature limit of the ther- 
mal conductivity is universal but depends on the slope 
parameter 



(49) 



The dependence of the universal value of the thermal 
conductivity on the slope of the excitation gap reflects 
the importance of the low-energy continuum states with 
e ^ Aq in the formation of the zero-energy impurity 
bound states. Note, that this result is valid for s-wave 
scattering and 7 ^ Aq with corrections typically of the 
order (7/Ao)^. A similar calculation for the in-plane con- 
ductivity gives Re(T||(T, w ^ 0) ~ 2e^Nf /(it /jAq), in 
agreement with Lee for the standard rf-wave model with 
= 2. 

For the three pairing states listed in Table I containing 
a line node in the basal plane (polar, hybrid-I, hybrid-II). 
the relevant angular average for heat flow in the basal 
plane is {j^Vf^ (A^ -|- 7^)"^^^) with Vf = Vf Vf , which be- 
comes i • d-d sin^ 1^ j^(A'^ + . The integral is 
dominated by the contribution near the line node, i.e., 
i9 7r/2. For 7 ^ Aq, we linearize the gap in the neigh- 
borhood of the line node, A(i9) ~ //Ao(7r/2 — and 
the integral reduces to rf??' 7^ [(//Ao«?')^ + 7^]"^''^) 
again leading to a universal result for 
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TABLE I. Symmetry classes, order parameters, and the asymptotic values of the thermal conductivity tensor. Note, that 
we have neglected the vertex corrections only for kj_ of the polar state. 



Pairing State Symmetry A(pj) Nodes ''ll(^) [ —m^^ ] '^-l(^) [ —m^^ 

T \ 3 ^ I T \ 3 ^ 

Class (Group) ^ ^ ^ 





Big (D4h) 




4 linear line nodes 


2 


— 




polar 


Alu (Deh) 


^PJ. 


1 linear line node 


1 


(tTa^) 


yAp 


liybricl-I 




Plz (PJx + ^Pjy ) 


2 linear point nodes 
+ 1 linear line node 


1 


~' 




hybrid-II 


E2u (Deh) 


^PtJPt:c + iPtyf 


2 quadratic point nodes 
-|- X linear line node 


1 


1 




2mAo 


2m2Ao 




hybrid-III-A 


Eig (Dgh) 


PJziPj^ + VtyKPtI + Pty) 


2 cubic point nodes 
+ 1 linear line node 


1 

2mAo 


0.47 / M3A0 ^ 

M3A0 V y ) 


1 


hybrid-III-B 


B2g+ «Big (Dgh) 


PtJPt:c + iPtyf 


2 cubic point nodes 
+ 1 linear line node 


1 

2mAo 


0.47 / M.iAn ^ 
M3A0 V T ) 


1 




Biu (Dgh) 


^Im(p/^ + ipjy'f 


2 cubic point nodes 
+ 3 linear line nodes 


3 

2mAo 


~10 / M3A0 ^ 
M3A0 V T ) 


1 

1^ 



TT" o 1 



«:||(T)~-#;t;/T 



(50) 



For heat flow along the c-axis the results for these same 
pairing states differ signiflcantly. First consider the 
hybrid-I state. The angular average now reduces to 
\ ■ /J^(i«?sin«?cos^ «?7^(A^ + 7^)~^/^. The gap opens up 
linearly at the positions of the point nodes at the poles. 
In the limit 7^0, the integrand diverges as \d — 
near the line node, and as near the point nodes. Thus 
the integral is dominated by the contribution from the 
point nodes (note, however, that 7 is mainly determined 
by the line node). Linearizing the gap near the point 
node A(«?) ~ //iAo«?, we obtain for the low-temperature 
limit of the c-axis component of the thermal conductivity 



1 



(51) 



which is nonuniversal, and generally much less than K|| 
by a factor of order (7/A0). To differentiate gaps with 
point nodes of different orders we classify them by their 
flrst nonvanishing derivative at the nodal points [i.e., for 
an n-th order point node, A ~ /J-n^o ^"]- 

In the case of the hybrid-II gap the c-axis transport 
is more subtle. The gap opens quadratically with angle 
near the point nodes at the poles, i.e., A(«?) ~ fi'j^oi^^, 
and once again we obtain a universal result for the zero- 
temperature thermal conductivity. 



-Njv]T 



1 



2//2A0 



(52) 



The ground state for the E2t, model^^ of UPts is an ex- 
ample of a hybrid-II state. Thus, an important feature of 
this model is that both components of the thermal con- 
ductivity tensor have universal values in the limit T ^ 0. 
If we use the polynomial form for the E2t, order parame- 
ter, A(p}) ~ Pfz{Pfx + Wfy)^! then the slope of the gap 



at the line node and the curvature of the gap at the point 
node are identical, ji = = 3-^3/2. Thus, for a spher- 
ical Fermi surface we flnd that all the eigenvalues of *k 
are identical in the limit T ^ 0. This result is consistent 
with the result reported in Ref. 47. However, the isotropy 
of peculiarity of the polynomial basis functions 

for the F,2u representation, and of course the spherical 
Fermi surface. 

Finally, we consider a gap with cubic point nodes 
at the poles, i.e., A(«?) ~ /isAo-d^. An example is 
the hybrid-III-A state with Ei^ symmetry, A(pf) = 

16/(3^3) Ao Pf zip} +p)y ) (Pfx +iPfy ) , which has a linear 
line node in the basal plane, but also cubic point nodes at 
= 0, TT. The Big and states of a hexagonal crystal 
possess cubic point nodes along the c-axis if one assumes 
an analytic expansion of A in terms oipj . The odd-parity 
Bit, and B2t, states also possess cubic point nodes if we re- 
strict the spin quantization axis to c? = c; however, more 
general spin states do not possess cubic point nodes. It 
is also worth noting that cubic point nodes are expected 
for a large number of superpositions of two ID represen- 
tations, as in the "accidental degeneracy" models®"'^ for 
UPts [c/, Refs. 52 and 39]. There is an important differ- 
ence between the 2D hybrid-III-A state and the various 
ID representations of the hexagonal group. The cubic 
point nodes of the ID or mixed symmetry ground states 
are generally combined with line nodes connecting the 
point nodes at opposite poles [An exception is the de- 
generate B2J -|- iQig (hybrid-III-B) state]. Thus, there is 
a higher density of excitations for the ID states than for 
the hybrid-III-A state in the vicinity of the cubic point 
nodes. The large density of gapless excitations in the 
vicinity of a cubic point node leads to a relatively large 
nonuniversal value for the low-temperature limit of the 
thermal conductivity in the c direction. For the case of 
a pure cubic point node we obtain 
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-Njv]T 



27rV37r 



27//3Aor(|)r(|) 



7 



(53) 



in the limit T ^ 0. Note, that the slope of the c-axis 
thermal conductivity is enhanced relative to the univer- 
sal limit of the quadratic point node by (//sAq/y)^''^. Of 
course this value is limited at large 7 by the normal state 
value. 

Although we have so far assumed an isotropic or cylin- 
drical Fermi surface, many of the results are more gen- 
eral, or simply extended to include uniaxial anisotropy. 
Also, one can generalize the results to an arbitrary 
anisotropic Fermi surface with anisotropic impurity scat- 
tering of the form = ri(pf)ri(pj'), where ri(pf) is 
any basis function with the full symmetry of the Fermi 
surface. In addition, for the zero temperature limit only 
the density of states and Fermi velocities at the appro- 
priate nodes are involved. Thus, the results for the zero 
temperature slopes of the thermal conductivity are easily 
extended to include Fermi surface anisotropy by replac- 
ing NfVj by the value of this quantity at the position of 
the node, or the relevant one-dimensional average in the 
case of a line node. 



B. Low Temperature Corrections 

Impurity scattering in an unconventional superconduc- 
tor with line nodes leads to a finite density of zero-energy 
excitations. The bandwidth of these impurity bound 
states is of order 7. The leading order finite temperature 
corrections to the transport coefficients are obtained by 
a Sommerfeld expansion of Eqs. (30) and (31) for the 
thermal and electrical conductivities. The key point is 
that the impurity-renormalized excitation spectrum has 
according to the symmetry relation (32) a low-energy ex- 
pansion of the form 



~R,A 



(e) ±i(j + he^) + ae , 



(54) 



with real coefficients a,h and < 7 ^ Aq. Expanding 
the integrands in Eqs. (30) and (31) to 0[e^] gives the 
Sommerfeld expansion for the components of Kij and aij , 



Kii(T) 



Njvf 



dee^sech^ ( ^] < 7^/3/2 



4T2 ./ """" V2TJ 
'5 



+ € 



Re aii(T, uj^O) 



27&^3/2+ (2«'7'-3&7')4/2 
-^«V/7/2]|, (55) 

e^Njvf 



AT 



desedi"^ (2^) y^^l"^ 
•15 



27&^3/2+ (yaV -3&7')4/2 

(56) 



—a 7 77/2 



where ly = ^^^^(A^ -|- 7^) ''^ for the in-plane compo- 
nents and ly = ^t)J^(A2 -|- 7^)"*^^ for the c-axis compo- 
nents. The different numerical coefficients for the terms 
involving are due to the difference in the coherence 
factors for electrical and thermal transport. For = 3/2 
these integrals have been evaluated: 73/2 = 2/(7rAo7^) 
for the dx2-y2 state and 73/2 = l/(2Ao7^) for the states 
with a line of nodes at i9 = 7r/2. We find in all cases 
/5/2//3/2 = 272/3 and 77/2/73/2 = 87VI5. It is re- 
markable that h always drops out. After performing the 
e integrals — Jde e'^"sedi\e/2T) = &„7r2"T2"+i, with 
60 = 4, &i = 4/3, &2 = 28/15 — we obtain for the dj;2_y2 
state 



3 TTfiAo 



15 72 



and 



Re(T||(T,w^O) 



,2^ 

7r//Ao 



7r2 a'T^ 
3 72 



(57) 



(58) 



to leading order in aT /■j. For the other pairing states 
(polar, hybrid-I, hybrid-II) the in-plane transport coef- 
ficients are obtained from Eqs. (57)-(58) by multiplying 
by 7r/4. The coefficient a is strongly dependent on the 
phase shift. For resonant scattering a = 1/2, indepen- 
dent of the specific pairing state. Note, that Eq. (58) 
agrees with the result of Hirschfeld it et al. in the res- 
onant limit. In the Born limit a = 7r//AoTo/2 for the 
dj;2_y2 state, and a = 2//AoTo for the other states (ignor- 
ing the special case of quadratic or cubic point nodes). 
Since we assume AqTo ^ 1, we always have for weak 
scattering a ^ 1. 

Finally, note that the finite temperature correction to 
the Wiedemann-Franz ratio becomes 



L(T) 



T Re (Til (T, 



■0) 



Ls 1 



27r2 a^T- 
15 72 



(59) 



which increases with temperature for T ^ T* ~ 7. This 
behavior arises from two sources: (i) the density of states, 
which is finite at e = with A^(0) ~ Nf (7/A0), depends 
strongly on energy for e ^ 7, and (ii) the difference in the 
coherence factors for thermal and electrical conduction, 
which also depend on e. Note that if scattering is weak, 
or if the material is very clean, then the very low tem- 
perature regime may be difficult to achieve in practice. 



V. NUMERICAL RESULTS 

More detailed information can be obtained from nu- 
merical evaluations of the transport coefficients over the 
full temperature range below Tg. The numerical results 
reported here were obtained by computing the equilib- 
rium propagator, self-energy, and order parameter self- 
consistently for the four pairing models — dx2-y2, polar, 
hybrid-I and hybrid-II — then using these results as input 
to numerically evaluate Eqs. (30)-(31) for the transport 
coefficients. We assumed a spherical Fermi surface except 
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for the 2D dx2-y2 state for which we used a cylindrical 
Fermi surface. 

Our numerical results agree with those of previous 
authors^'*''*®''*^''*^'®^^®® in those cases where a direct com- 
parison is possible. 

Figure 1 shows the results of the in-plane thermal con- 
ductivity for the four pairing states in the resonant scat- 
tering limit, i.e., for a normalized scattering cross section 
of (T = sin^ (5o = 1. All curves exhibit the qualitative 
behavior of a superconductor with lines of nodes as dis- 
cussed by many authors. It is remarkable that the curves 
for the 2D dj;2_y2 state and the 3D E2t,-state are essen- 
tially identical. 
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FIG. 1. Thermal conductivity Ky vs. temperature for un- 
conventional superconductors in the unitarity limit {u = 1) 
with a dimensionless scattering rate a = 0.01. The different 
pairing states are: (i) d^2_y2, (ii) polar, (iii) hybrid-I, and 
(iv) hybrid-II as described in the text. For comparison the 
result for an isotropic BCS superconductor is shown. Panel 
(a) displays the in- plane or basal plane Ky normalized to its 
value at Tc- Panel (b) displays the ratio Ky/T normalized 
to its value at Tc- Note the finite intercept of Ky/T for the 
unconventional pairing states. 

In Fig. 2 we plot for a dj;2_y2 pairing state 
[k\\(T)Tc / K\\(Tc)T] vs. T/Tc for several (normalized) 
scattering cross sections a = sin^ 60 and (normalized) 
scattering rates a = I / (2'itTcoTo) ■ A consequence of the 
universal limit for T ^ is that the ratio 



K\\(T)T, 



1 



2 Tco 



T™ K|| (Te)T 7rroAo(0) Ao(0) 



(60) 



scales linearly with the scattering rate parameter a, pro- 
vided a is significantly less than the critical pair-breaking 
value (ttcr ~ 0.28), and is independent of the scattering 
strength, a. At T ^ T* and for weak scattering Arfi ei 
al.^^ have shown that the ratio [k||(T) Tc]/[k||(Tc) T] cx 
(I — a) strongly depends on the scattering phase shift 60- 
This explains (i) the sudden drop of [k||(T) 7"c]/['^||(7c)] T 
in Fig. 2(a) at ultra-low temperatures for weak scat- 
tering, where the universal limit is achieved only for 
T ^ T* ~ Aoexp(— 1/a), and (ii) the scaling of the zero 
temperature intercept in Fig. 2(b). To address the vari- 
ous power law behaviors of Ky (T) in different temperature 
regions and for different scattering rates and scattering 
cross sections, we show K||(T)/k||(Tc) in Figs. 2(c,d) in a 
log-log plot for the same parameters as in Figs. 2(a,b), re- 



spectively. The temperature dependence of the electronic 
thermal conductivity obeys a variation above a crit- 
ical temperature T* ~ 7 in clean superconductors and 
in the strong scattering regime. Below T* it approaches 
the limiting T behavior. Weak scattering leads to an ap- 
proximately linear temperature dependence over a large 
portion of the temperature range. However, the ratio 
K||/T changes drastically in clean superconductors below 
the exponentially small crossover temperature T* , where 
it approaches its linear low-temperature asymptote. 
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FIG. 2. Thermal conductivity of a 2D rf^2_j,2-wave super- 
conductor. The temperature dependence oi n\\yF)Tcl n\\{Tc)T 
vs. T/Tc is displayed in (a) for different scattering cross sec- 
tions, (J = sin^ ^0, and for a fixed scattering rate a = 0.01. In 
panel (b): the temperature dependence is plotted for differ- 
ent scattering rates a in the unitarity limit (a = 1.0). Panels 
(c) and (d) are the same calculations as panels (a) and (b), 
respectively, plotted on a log-log scale in order to exhibit the 
low temperature power laws and crossovers. The thick dotted 
lines are guidelines to T and T'^ power laws. 



As T ^ 0, we indeed find the universal behavior as 
discussed in the previous section. To show the approach 
to the universal limits at low temperatures, we computed 
the electrical and thermal conductivity at low and ultra- 
low temperatures for an intermediate scattering rate a = 
0.1, chosen because the temperature range of universality 
is exponentially small for weak scattering cx exp(— 1/a). 
The results are shown in Fig. 3 where the electrical and 
thermal conductivities have been normalized by their cor- 
responding universal limits, (Tq = (T||(T ^ 0, w ^ 0) and 

Ko(T) = 7""^(0). The corresponding Lorenz ratio is 
shown in Fig. 3(c). Notice the logarithmic scale in tem- 
perature. In agreement with the analytical results, nu- 
merical calculations show that Re(T||(T, w^O), K||(T)/T 
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and L(T) increase with temperature near T = 0. The 
crossover temperature to the universal regime is expo- 
nentially small for small phase shifts, and the Lorenz 
number quickly drops for these smaller phase shifts after 
the initial rise. 
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FIG. 3. Transport coefficients Re (7|| (T, a; ^ 0) and K||(T) 
vs. T of a 2D rf^2_j,2-wave superconductor for different scat- 
tering cross sections (a = ... 1) and a fixed intermediate 
(normalized) scattering rate a = 0.1. (a) Tlie normalized elec- 
trical conductivity Re (7||(T, a;)/(7o at a very small frequency 
ui = 10~^Ao(0). (b) Tlie normalized tliermal conductivity 
K|| (T)/fio (T). (c) Tlie normalized Lorenz ratio i(T)/is. Note 
tlie significant difference in L{T)/Ls for unitarity vs. Born 
scattering. Inset: Blow-up of the universal behavior at ul- 
tra-low temperatures showing the approach to is at T ^ 0. 




FIG. 4. The normalized Lorenz ratio L{T)/Ls of a 2D 
rf^2_j,2-wave superconductor for different (normalized) scat- 
tering rates a in the Born (a = 0) and unitarity (a = 1) 
limit. L{T)/Ls clearly separates between scattering in the 
weak limit (lower four curves) and the strong limit (upper 
four curves). 



In Fig. 4 we show the Lorenz ratio, L(T)/Ls, over the 
full temperature range for scattering in both the weak 
(Born) limit and the resonant (unitarity) limit. The de- 
viations from the Sommerfeld value, Ls , are clearly sepa- 
rated between the Born and unitarity limits. This effect 
is most pronounced in nearly pure systems. At temper- 
atures close to Tc, and in the clean limit with strong 



scattering, the Lorenz ratio is slightly reduced due to a 
small coherence peak in the electrical conductivity. In 
the Born limit L(T)/Ls is always less than the normal 
state limit, which is opposite to that for unitarity scatter- 
ing. At sufficiently low temperatures all curves converge 
to the same universal limit. 

In superconductors with a large concentration of (res- 
onant) scatterers the Wiedemann-Franz law is fulfilled 
throughout the entire temperature range, a result which 
is obvious from the fact that nonmagnetic impurities lead 
to pair-breaking in unconventional superconductors, and 
as the impurity lifetime tq approaches the critical pair- 
breaking value, the transport properties approach those 
of the normal metal. 



VI. COMPARISON WITH EXPERIMENTS 

A. High Tc cuprates 

Several measurements of the thermal®^'®® and 
electrical®^'®*^ conductivity on single crystals of high Tc 
cuprates have been reported. In-plane thermal con- 
ductivity measurements show the presence of a 'low- 
temperature' term, which has been interpreted as 
boundary scattering of phonons on crystal faces, as well 
as a linear term at very low temperatures, which has been 
attributed to uncondensed charge carriers [for a review 
see Ref. 1]. 

In Table II we list the material parameters for the 
cuprate superconductors which we used to estimate the 
slope of the thermal conductivity at zero temperature. 
In the analysis we used the Drude plasma frequency to 
determine NjvJ , i.e., oj'^ = Ai:e'^ NjvJ , the weak-coupling 
(i-wave gap ratio, Ao(0) = 2.14A;bTc, and the universal 
lifetime ta — l/(7rAo) to rewrite the universal limit of 
the electrical conductivity as 



47r 



(61) 



This enables us to estimate the slope of the thermal con- 
ductivity from the WF law. 



lim (K(T)m 



(62) 



For Y-Ba-Cu-Oe 95 the conductivity at T ^ esti- 
mated from microwave experiments is Ua ~ 0.5 (//SI m)~^ 
along the a-axis and (Tj k, 0.7(//Slm)~^ along the b- 
axis.®"'^ These values are close to the universal value 
obtained from Eq. (61) and Table II, (Tq ~ 0.5 — 
0.6 {jiVl m)~^, and provide reasonable agreement between 
the slope of the thermal conductivity obtained from the 
Wiedemann-Franz ratio (62) and the experimental val- 
ues listed in Table III.®^ However, the experimental coef- 
ficients vary by as much as a factor of ~ 8 among differ- 
ent samples of Y-Ba-Cu-Oe 95. The theoretical values for 
limy^o ^^w/T lie in the range of experimental values ex- 
cept for La-Sr-Cu-0. Thus, if La-Sr-Cu-0 is a rf-wave su- 
perconductor resolution of this discrepancy would require 
weak-scattering. In which case it would be very difficult 
to determine the universal zero temperature slope of Ky, 
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because the region where Ky shows universal behavior is 
exponentially small. Thus, for Born scattering the low 
temperature extrapolation of Ky/T will overestimate the 
universal limit. 



TABLE II. Material 
cuprate samples. 



parameters of optimally doped 



Compound 


Tc (K) 




(eV) 


Refs. 


La-Sr-Cu-O 


35 


0.8 


- 0.9 


68, 69 


Y-Ba-Cu-O 


92 


1.4 


- 1.5 


61, 70, 71 


Bi-Sr-Ca-Cu-O 


86 


1.1 


- 1.2 


70, 72 



TABLE III. Slopes of the low-temperature thermal con- 
ductivity for the cuprates. The theoretical value is the uni- 
versal slope calculated for a d^2_y2 pairing state in the limit 
of resonant scattering (see Sect. VI). 





La-Sr- 


■Cu-O 


Y-Ba-Cu-O Bi-Sr-Ca-Cu-O 


(mW/K^ m) 








theor. 


10 - 


■ 13 


11-14 8-10 


exp. 


100 - 


■ 150 


5 - 12, 40 ~ 6 


Refs. 


73, 


74 


1, 58 1 



B. The heavy fermion superconductor UPts 

For UPts and heat flow in the basal plane the ther- 
mal conductivity in the normal state obeys k„(T)/T = 
l/{a + bT^), with a = 0.25mK2W-i and b = 
1.0 mW~^, and it is believed that the thermal con- 
ductivity is almost entirely electronic,^® and that the 
h term arises from electron-electron scattering. If we 
take the experimental data reported by Lussier, Ell- 
man and Taillefer^® for UPts with Tg ~ 0.5 K, an av- 
erage Fermi velocity, Vf = 5.5km/s, and a transport 
mean-free path of /q = 220 nm [i.e., a transport time of 
To = 40 ps] , and combine it with the theoretical value for 
limT-^ K||(T)/T 0.6 ^^^(o) K„(Tc)/Tc, we can predict 
the universal low-temperature slope. For the F,2u state 
we estimate aslope of limy.^ K||(T)/T 0.1 k„(Tc)/Tc 
2mWK~^cm~^, while the corresponding ratio for the 
Eig state is smaller by 20%. These estimated theoreti- 
cal values of Ky/T at T = 0, are of the same magnitude 
as the experimental data at T 0.1 ~ 50 mK with 
Ky/T|e^p ~ 2mWK-2cm-i.29 

To distinguish these different pairing states it is neces- 
sary to probe a symmetry dependent quantity as the non- 
linear Meifiner effect, ^^'^"^ or the -l-phase shift Joseph- 
son effect. Here we analyze the anisotropy ratio of 
the thermal conductivity. The experimentally observed 
anisotropy ratio 



lim — — 

T^O (ki/k||)„ 



(63) 



approaches r-a — 0.43 for T 0.^® For an ellip- 
soidal Fermi surface the normal state anisotropy, e.g., 
Vf±_ Vf\\, drops out of the anisotropy ratio in Eq. (63). 
Thus, we can use the formulas derived for an isotropic 
Fermi surface; the anisotropy is then determined by the 



anisotropy of the superconducting gap as deflned by the 
gap anisotropy parameters, //i, 1^2, which measure the 
slope or curvature of the gap at a line or point node on 
the Fermi surface.®'* 

The anisotropy ratios for the various order parame- 
ters are obtained from Table I. The Ei^ and F,2u models 
have been discussed by many authors as models for the 
low temperature phases of UPt3.®5,46, 66,67,39,11 the 
Eij (hybrid-I) state the anisotropy ratio is nonuniversal 
and given by 



27/" 



(Eig) . 



(64) 



For the standard Ei^ order parameter, A(p) = 
2AoPfz(Pfx + ipjy), we have jj, = ni = 2, and hence Va — 
7/A0. Assuming the unitarity limit and = 0.1 T^, 
we obtain Va ~ 0.2. We can account for the experi- 
mental value of Va by adjusting jj, and/or jii. Although 
Ft, = 0.1 Tc is consistent with Ref. 28 and other nor- 
mal state measurements, it predicts a flattening of kj 
for T ^ 0.2 Tc, which is experimentally not observed. A 
smaller scattering rate. Ft, = 0.01 Tc, does a better job in 
accounting for the low temperature behavior of K|| , hence 
we obtain Va ~ 0.04. Thus, it is not possible to account 
for the experimental anisotropy, Va — 0.43,^® without a 
very large ratio fi/fif, which seems unphysical. 

However, for the F,2u (hybrid-II) state the anisotropy 
ratio becomes 



ra ~ /U//U2 (E2u) : 



(65) 



which is universal and independent of 7. Thus, we can 
easily flt the experimental value of fa for T ^ by choos- 
ing /u//U2 ~ 0.43, and then adjusting 7 to obtain the best 
flt to the low temperature behavior for K||(T). 

Finally, we note that for the cubic point nodes the 
theoretical prediction for the normalized anisotropy ra- 
tio would be 



1^3 



(hybrid - III) 



(66) 



Therefore we generally expect fa > 1, which clearly dis- 
agrees with experiment unless one chooses a very large 
IJ,s. Barring this unlikely scenario we can conclude that 
the anisotropy ratio rules out cubic (or higher order) 
point nodes for the gap in UPts. 



VII. CONCLUSION 

In conclusion we have examined in detail the low- 
temperature behavior of the thermal and electrical con- 
ductivities and the Lorenz ratio for superconductors with 
line nodes. We considered several unconventional pair- 
ing states and found for some of the eigenvalues of the 
thermal conductivity tensor a universal value as T^O, 
similar to the electrical conductivity. Furthermore, we 
showed that the Wiedemann-Franz law is restored below 
a crossover temperature T* . In clean systems we have 
T* ^ Tc, while in 'dirty' superconductors the crossover 
temperature can be an appreciable fraction of Tg. 
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The estimates of the universal slope limT_^ k/T, which 
we derived for the cuprates based on a dj;2_y2 order pa- 
rameter, and for UPts, with a hybrid-I or -II gap, are 
comparable with experimentally reported values, except 
for the La-Sr-Cu-0 materials. Assuming that La-Sr- 
Cu-0 is an unconventional superconductor, the differ- 
ence between theoretically and experimentally extrap- 
olated values might be attributed to weak scatterers 
in this system. In UPts the anisotropy ratio, ra = 
limT^o('^J_/'^||)/('^J_/'^||)n is consistent with an F,2u gap 
and a universal value for ra, or with an Ei^ gap and 
a nonuniversal value for ra determined by the impurity 
concentration. Further experiments on UPts with con- 
trolled impurity concentrations should easily distinguish 
these two models 
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APPENDIX: SOLUTIONS TO THE LINEARIZED QUASICLASSICAL TRANSPORT EQUATIONS 



The deviations from local equilibrium, b'g 
linearized equations 



X _ -x -x ^^x ^ _ ^-^j^ X ^ {R,A,K], satisfy the 



9 - 9b ) anc 



69 



R,A i,R,A 



and 



'R „ 



Wffo + \9o , '^ext + OCT ' 



idg^+ {aext + 6a^) o g^ + o g^ 

-9o ° i^^ext + l5(T^) - ff^ O ,5(7^ , 



(Al) 



(A2) 



where h^'"^ = efs — ir^'^ and d = Vf ■ V . Note, that the equation for the Keldysh propagator, 6g^ , is coupled to 
the deviations of the retarded and advanced functions, b 'g^'^. We decouple these equations by using the equilibrium 
relations, Eqs. (23) and (24), and by introducing Eliashberg's anomalous propagator, bg'^ , and self-energy, ba'^ , defined 

by 



8 'g'^ = 8 'g^ o $0 - $0 o 8 'g^ + bg"" 
o $0 - $0 o + 6a"- 



(A3) 
(A4) 



After eliminating idg^ and idgQ using Eq. (Al), the transport equation for 6g^ is transformed into a transport 
equation for Sg", 

h^ob'g'' - bg"" oh^ 
+ g§ o (i5$o) - (i5$o) o + o 6a" - 6a" o g^ 

+ [(Text , ^o]o °9o +9o ° ['^0 , ^ext]o 



(h^ + a^) o <!>o - <!>o o (h^ + a^) 



• 6r 



6g' 



The terms on the right hand side of Eq. (A5) vanish identically because 

tA,-A_'iR.'^R_ 

n + (Tg — /i + (Tg — era . 
Thus, the equation for the anomalous propagator, 6g" , becomes 



(A5) 



(A6) 



h^o6g" - 6g" o = 
(i5$o) o9o-9oO (i5$o) + 6a" o g^ - gR o 6a" 

- [cText , $o]o °9o-9o° [*0 , (rext]o ■ (A7) 

The transport equations for 6'g^'^'" are solved by noting that the local equilibrium propagators have the form 

R,A 

(A8) 



^R,A 
9o 



i^'^T3 - A 



CR.A 



and that 



where 



j^R,A ^ ^jR,A ^R.A ^ ^R,A | 



C 



R,A 



1 



(A9) 



(AlO) 



and the functions e^'^, A^'^, and D^'^ are defined in terms of the equilibrium self-energy 

' R,A 



(Tg — (e - £ ' )r3 -I- A 



(All) 
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Finally, Eq. (A7) is solved by using the normalization condition for Sg", 

9^06^+6^09^=0, (A12) 

and the relations for /i^'^ in Eq. (A9), to move Sg" to the right in the second term on the left side of Eq. (A7). Then 
(A7) turns into a form which can be solved by matrix inversion, 

6g-= [C^g^ + D-_]-' o (^{id<fo) o - o {id<fo) 
+ 6(7" og^ - g^ o 6a"- 
- [&ext , $o]o og^ -g^o [$o , -Te^tlo) > (A13) 

where 

Cl(pf;e,L,) = C''(pf-,e+) + C^(pf-,e.), (A14) 
{pj ■e,oj) = D^ipj ■,e+)-D^ {pj ; e_ ) (A15) 

with e± = e ± w/2. Using the normalization condition (26) we write the inverse matrix as 



t^.r + i^J =-7l^^- (A16) 

Furthermore it is useful to introduce 

Cf^{pj;e,uj) = C^'^{pj;e+) + C^'\pj;e.) , (A17) 
= D^'^{pr,e+) - D^'^ipj-e.) , (A18) 

and to obtain the solutions for 6'g^'^ by analogous steps, 

8g^'^ = [Cf^g^^^ + D^J^\ o + + 8a^'^ , ,f '^J J . (A19) 

When combined with b'g'^ in Eq. (A13) according to (23) we obtain the general result for 6'g^ , 
Sg"" = [C^go +D^r'o (^-idg^ + [&e.t + Sa"" , g^] „^ o $0 
-$0 o [C^g^ + D^]-\ (^-idg^ + [a,., + 6a^ , g^] }j 
+ [Clg^ + Dl]-' o (^iid<fo) og^-g^o (i5<I>o) + 6a" o g^ - g^ o 8a" - 

{(Text , $o]o og^ -g^O [$o , <Te^t]„^ . (A20) 



-1 



Cg^-^ - D 
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